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ABSTRACT 

We present observations of the recently discovered supernova 2008iz in M82 with the VLBI High Sensitivity Array at 22 GHz, the 
Very Large Array at frequencies of 1.4, 4.8, 8.4, 22 and 43 GHz, and the Chandra X-ray observatory. The supernova was clearly 
detected on two VLBI images, separated by 11 months. The source shows a ring-like morphology and expands with a velocity of ~ 
23000 km s"'. The most likely explosion date is in mid February 2008. The measured expansion speed is a factor of ~2 higher than 
expected under the assumption that synchrotron self-absorption dominates the light curve at the peak, indicating that this absorption 
mechanism may not be important for the radio emission. We find no evidence for an asymmetric explosion. The VLA spectrum shows 
a broken power law, indicating that the source was still optically thick at 1.4 GHz in April 2009. Finally, we report upper limits on the 
X-ray emission from SN 2008iz and a second radio transient recently discovered by MERLIN observations. 

Key words. (Stars:) supernovae: general, (Stars:) supernovae: individual: SN2008iz, Radio continuum: general. Galaxies: individual: 
M82 



1. Introduction 

• Radio supernovae (RSNe) are rare events and difficult t o study. 
I So far only about two dozen have been detected (e.g Weile r et aH 
. I2OOI and most of them are quite distant and rather weak. 
' Only few radio supernovae have been imaged with Very Long 
Baseline Interferometry (VLBI) techniques and for only four has 
it been possible to study the evolution of the expanding shell 
(SN 1979C, SN 1986J, SN 1987A, SN 1993J). The best studied 
. radio supernova so far is SN 1993 J in M81. Following the ex- 
■ pansion of the supernova allowed many different phenomena to 
; be studied dM arcaide et al. 1997, 2009a; Bietenholz et al. 2001, 
'. I2003t iPerez^ on-es et al.l i200lL |2002a; Bar tel et al.ll2002l l2007l) . 
' including a measurement of the expansion speed, the decelera- 
' tion of the shock front, and the proper motion of the supernova 
shell, for which a limit was obtained. 

The recent di scovery of a new brig ht radio supernova in 
M82, SN 2008iz, (iBrunthaler et al.l2009allbh at a similar distance 
as SN 1993 J offers the rare opportunity to study the evolution of 
another supernova in great detail and to make a comparison to 
SN 1993J. So far, SN 2008i z was only detected in t he radio band, 
with the VLA at 22 GHz ( Brunthaler et"ani2009ah . MERLIN at 
5 GHz dg eswick et al. 2009), and the Urumqi telescope at 5 GHz 
(Marchili et al. 2010). There are no reported detections in visi- 
ble light and Fraser et al. (2009) report only a non detection in 
the near infrared on 2009 June 1 1 . The non-detections at other 
wavebands indicate that the supernova exploded behind a large 



gas or dust cloud in the central part of M82. Thus, it has not been 
possible to classify this supernova. However, since type la super- 
novae are not known to show strong radio emission, SN 2008iz 
is most likely a core collapse supernova, i.e. type Ib/c or II. 

Here we present the first VLBI images of SN 2008iz taken 
~2.5 and ~13.5 month after the explosion, a radio spectrum (at 
an age of ~14.5 month) from 1.4 to 43 GHz and Chandra X- 
ray observations. Throughout the paper, we adopt a distance of 
3.6 Mpc (based on a Cepheid distance to M81 determined by 
Freedrnanetal ■111994 . 

2. Observations and data reduction 

2.1. High Sensitivity Array observations at 22 GHz 

M82 was observed using the High Sensitivity Array includ- 
ing the NRACQ Very Long Baseline Array (VLBA), the Very 
Large Array (VLA), the Green Bank Telescope (GBT), and the 
Effelsberg 100m telescope under project code BB255 on 2008 
May 03 and 2009 April 08. The total observing time at each 
epoch was 12 hours. We used M81*, the nuclear radio source 
in M81, as phase calibrator and switched between M81*, M82, 
and 3 extragalactic background quasars every 50 seconds in the 
cycle M8 1 * - 0945+6924 - M8 1 * - 0948-1-6848 - M8 1 * - M82 

' The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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- 1004+6936-M81*, yielding an integration time of ~100 min- 
utes on M82. J1048+7143 was observed every ~30 minutes to 
phase-up the VLA. Furthermore, DA 193 was observed as fringe 
finder. The data were recorded with four 8 MHz frequency bands 
in dual circular polarization, with Nyquist sampling and 2 bits 
per sample (i.e., a total recording rate of 256 Mbit s"'). 

Before, in the middle, and at the end of the phase referencing 
observations, we included geodetic blocks, where we observed 
18-21 bright sources from the International Celestial Reference 
Frame (ICRF) at 22 GHz for ~75 minutes to measure the tro- 
pospheric zen ith delay offsets at each antenna (for a detaile d 
discussion see Reid & Brunthaler 2004: Brunthaler et al .1 120051) . 
The geodetic blocks were recorded with 8 IFs of 8 MHz in left 
circular polarization. The first IF was centered at a frequency of 
22.01049 GHz, and the other seven IFs were separated by 12.5, 
37.5, 100.0, 262.5, 312.5, 412.5, and 437.5 MHz, respectively. 
Here, only a single VLA antenna was used, which could observe 
only five IFs at 0.0, 12.5, 37.5, 412.5, and 437.5 MHz relative to 
22.01049 GHz. 

The data were correlated at the VLBA Array Operations 
Center in Socorro, New Mexico. In the observation on 2008 May 
03, the data were correlated at the position of known water maser 
emission (09*'55"5P.38702, -)-69°40'44".4676, J2000) with 128 
spectral channels per IF and an integration time of 1 second. 
The observation on 2009 April 08 was correlated at the position 
09*'55™5P.5500, +69°40'45".792, (J2000), close to SN2008iz 
with 16 spectral channels and an integration time of 1 second. 

The data reduction was performed with the NRAO 
Astronomical Image Processing System (AIPS) and involved 
standard steps. First we shifted the position of M82 to the posi- 
tion of SN 2008iz using CLCOR in the observation on 2008 May 
03. We applied the latest values of the Earth's orientation param- 
eters and performed zenith delay corrections based on the results 
of the geodetic-like observations. Total electron content maps 
of the ionosphere were used to correct for phase changes by 
the ionosphere. A-priori amplitude calibration was applied us- 
ing system temperature measurements and standard gain curves. 
We performed a "manual phase-calibration" using the data from 
J 1048+7 143 to remove instrumental phase offsets among the fre- 
quency bands. Then, we fringe fitted the data from M81* and 
performed an amplitude self-calibration on M81*. In the obser- 
vation on 2008 May 03, we also performed, first, a phase-self 
calibration and later an amplitude self-calibration on SN 2008iz. 
In the observation on 2009 April 08, the inclusion of the phased 
VLA increased the noise in the images significantly, possibly 
due to a problem with the phasing of the array (the VLA was 
in B -configuration). Thus, we flagged the VLA data from this 
epoch. 

2.2. VLA observation on 2009 April 27 (1.4-43 GHz) 

We observed M82 with the VLA on 2009 April 27 at 1.4, 4.8, 
8.4, 22, and 43 GHz. The total observing time was 4 hours. We 
observed with two frequency bands of 50 MHz, each in dual cir- 
cular polarization. 3C48 and J1048+7143 were used as primary 
flux density and phase calibrators, respectively. At 1.4, 4.8, and 
8.4 GHz, we used a switching cycle of six minutes, spending 
one minute on the phase calibrator and five minutes on M82. We 
repeated these cycles 5 times over the observations, yielding an 
integration time of ~25 minutes at each frequency. At 22 and 43 
GHz, we used a switching cycle of three minutes, spending one 
minute on the phase calibrator and two minutes on M82. These 
cycles were repeated 10 times during the observation, yielding 
an integration time of ~20 minutes at both frequencies. 
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Fig. 1. HSA contour plots of SN 2008iz at 22 GHz in May 2008 
(top) and April 2009 (bottom). Contours start at 0.3 mJy (top), 
0. 1 mJy (bottom) and increase with V2. The images are restored 
with a circular beam of FWHM 0.2 mas (top) and 0.5 mas (bot- 
tom). The original beam sizes were 0.31x0.23 mas at a position 
angle of -34° at the first epoch and 0.39x0.27 at a position angle 
of -40° at the second epoch. 



The data reduction was performed in AIPS and involved am- 
plitude calibration of 3C48 using source models. Then we cal- 
ibrated the phases using J 1048 +7 143 and made one phase and 
amplitude self-calibration on J 1048 +7 143. The calibration was 
then transferred to the target source M82. We performed one 
phase self-calibration on M82 at all frequencies except 43 GHz, 
where the source is too weak. 

2.3. Cliandra and Swift XRT observations 

The Chandra X-ray observatory (CXO) observed M82 on 2008 
October 4, and 2009 April 17 and 29. Each observation was 
taken with an exposure of about 18 ks. In these observations. 
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Table 1. Flux densities of SN 2008iz and two other supernova 
remnants in M82 from the VLA observation on 2009 April 27. 



SN 2008iz on 2009 April 27 



Source 


Fl.4 GHz 


F4.8 GHz 


F8.4 GHz 


F22 GHz 


F43 GHz 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


SN2008iz 


55.3±6.0 


29.5±4.5 


13.7±2.5 


6.2+0.7 


2.5+0.4 


44.01+596 




14.7±5.3 


15.0±2.0 


8.7+0.6 


3.6+0.4 


45.17+612 




8.5±1.6 


5.8+0.7 


3.0+0.9 


0.8+0.3 



3. Results 

3. 1 . VLBI images at 22 GHz 

We imaged the data of SN 2008iz from 2008 May 03 with a 
slightly super-resolved beam of 0.2 x 0.2 mas. The peak flux 
density was 5.9 mJy beam"', the total emission was 32 mJy and 
we achieved an image rms of 79 fiiy beam"'. The data from 
2009 April 08 were imaged with natural weighting and a circu- 
lar beam of 0.5 x 0.5 mas to have more sensitivity for extended 
emission. Here, the peak flux density was 0.36 mJy beam"', the 
total emission was 4.3 mJy, and we achieved an image rms of 41 
fijy beam"'. 

The supernova was clearly detected in both epochs (see 
Fig. [1]). The source is already clearly resolved in the observa- 
tion on 2008 May 03 and shows a ring like structure, typical 
for a radio supernova. In the following eleven months the source 
expanded and faded significantly. 

3.2. VLA radio spectrum 

M82 was imaged at all frequencies using only data from base- 
lines larger than 30 kA. This ensures that most of the extended 
emission in M82 is resolved out. Flux densities were extracted 
by fitting two-dimensional Gaussians to the images. 

The measured flux densities of SN 2008iz and two other su- 
pernova remnants that could be easily separated from the diffuse 
background emission are listed in Table [T] The errors are esti- 
mated by adding in quadrature the formal error from the fit to 
the images, the difference between peak and integrated flux den- 
sities, and an 5% error allowing for an error in the overall flux 
density scale. Furthermore, we added an additional 5% error at 
1.4, 4.8, and 8.4 GHz since we have more confusion from the 
extended emission at these frequencies. Note that the flux den- 
sity at 4.8 GHz is in good ag reement with the flux densities re- 
ported bv lBeswick et al.l (I2009I) . obtained with MERLIN a few 
days later (28.5± 2 mJy). 



the target is off" the optical axis by more than 3 arcmin, where 
the point spread function looks like an extended ellipse covering 
multiple pixels. 

M 82 was also observ ed with the Swift X-ray Telescope 
(XRT;|Gehrelsetal. 2004; Burrows et alJ2005h on 2007 January 
26, 2008 May 1, and 2009 April 25 with exposures of 4.6 ks, 5.0 
ks, and 4.7 ks, respectively. Due to the low resolution of the tele- 
scope, the collection of known X-ray sources looks point-like in 
the XRT data, so count rates from the conglomerate of sources 
were measured. The data suffer from heavy pile-up, which is 
accounted for by extracting count rates from an annular region 
around the central, piled-up location 0. 




10 
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Fig. 2. Spectrum of SN2008iz taken with the VLA on 2009 
April 27. Also shown are a single power-law fit (dotted line) 
and two broken power-law fits with a spectral index of 2.5 (solid 
line) and 4.5 (dash-dotted fine) in the optically thick part. 



The spectrum of SN 2008iz is shown in Fig.|2] First, we fitted 
a single power-law spectrum to the data. This gives a spectral 
index of -0.88+0.07. However, the fit has a large reduced 
value of 2.6, since the 1.4 GHz value is too low. Thus, we fitted 
the spectra also with a broken power-law. 



5(v) = 5„ - 



1 -e ^""Z 



(1) 



^ Given the pile-up and uncertainty in the XRT astrometry on the 
order of a few arcseconds, attempting to separate any supernova flux 
from the other flux is practically impossible 



where a and 6 are the spectral indices of the optically thin and 
thick parts of the spectrum. Sq and vq represent the maximum 
flux density and the peak frequency of the fitted spectrum. Since 
one data point in the optically thick part of the spectrum is not 
enough to fit the spectral index there, we made two fits, one with 
a value for a synchrotron self-absorbed spectrum {S =2.5), and 
a steeper free-free absorbed spectrum (5=4.5). The reduced 
values in both cases are now 0.9. The fit (using S = 4.5) gives 
a spectral index in the optically thin part of a = - 1 .08 ± 0.08 
and a turnover frequency of vq = 1.51 + 0.09 GHz. While the 
spectral index a is not affected by the choice of 6, vq changes 
slightly to 1.55 for S - 2.5. This indicates that the source was 
still optically thick, and brightening at the lowest frequencies in 
April 2009 (for comparision: SN 1 993J reached it's pe ak at 1.4 
GHz -500 days after the explosion: IWeiler et al.ll2002l) . 

3.3. X-ray upper limits 

The off-axis configuration and consequently the degraded angu- 
lar resolution, as well as the diffuse background, have largely 
decreased the sensitivity of the Chandra measurements. The de- 
tection limit of these observations is estimated from the total 
emission around the radio position in a region with the same 
size as the point spread function. SN 2008iz is located close to 
several variable ultraluminous X-ray sources. Ho wever, no emis- 
sion w as de tected at the position o f the supernova. iMuxlow et al.l 
(I2009h and lMuxlow' et al.l (l2010t) report the discovery of a sec- 
ond radio transient in M82 with the MERLIN telescope. This 
source appeared between 2009 April 24 and 2009 May 5 and is 
located at a position with a diffuse emission background. It is 
surrounded by a few point-like sources. There is no enhanced X- 
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ray emission at the location of this second radio transient, neither 
on 2009 April 17, nor on 2009 April 29. 

To compare our limits with the emission from SN 1993 J, 
we assume spectral properties that are similar to the ones of 
SN 1993J at a similar age, i.e., a thermal bremsstrahlung spec- 
trum with a te mperature of 1.05 keV, abundances from Table 
2, Column 4 of ^Zimmerma nn & AschenbachI (l2003b . and an ab- 
sorption column density of 5.4x10^^ cm"^ (see Sect. 15. lb . Then 
we get an 3 cr upper limit of 1.5x10^^ erg s"' in the energy 
range 0.3-2.4 keV. This is consis tent with the X-ray luminos- 
ity o f SN 1993 J at a similar age dZimmermann & AschenbachI 
12003 ) in the same energy range. 

The 3 cr sensitivity at the location of the MERLIN transient, 
assuming a column density of 10^^ cm"^ and a photon power- 
law index of 1 .7 (an approximation to a thermal bremsstrahlung 
spectrum with a temperature of 10 keV, as seen in SN1995N; 
iFox etal .'20Q0). is found to be about 1.2x10^** erg s^' in 0.3-2.4 
keV, which we take as the upper limit of the X-ray luminosity. 

The resultant count rates in the Swift XRT data (0.5-8.0 keV) 
are 0.60+0.02, 0.59+0.02, and 0.68+0.02 s"' on 2007 January 
26, 2008 May 1, and 2009 April 25, respectively. There is no 
clear rise in total flux from the pre-SNe observation to the May 
2008 post-explosion epoch at ~75 days, and while a ~15% in- 
crease in count rate is seen in the final observation, this is consis- 
tent with BeppoSAX observations of the central region of M82 
showing var iations on the ord er of 15-30% (2-10 keV) on hour 
time-scales dCaPDi etalJlf999h . Indeed, such intrinsic variations 
are also seen over the course of each epoch in the XRT obser- 
vations as well. Given the complications with the data and the 
intrinsic variability of the sources, strong constraints are diffi- 
cult. However, since no increase in flux larger than the intrinsic 
variability is seen, we can place an upper limit on the SNe X-ray 
flux of approximately 1.5x10'*' erg s"' (0.5-2.4 keV) on 2008 
May 1 (assuming the same model as for the Chandra observa- 
tions above). 



4. Expansion speed and explosion center 

4.1. Estimates of size and position 

Estimating a physical size from the VLBI data is not straight- 
forward. Different methods have been applied to d etermine the 
expansion curves of several radio supernovae (e.g. iBartel et"an 
[2002; Marcaide et al. 2009a.b). Estimating the source size in the 
image plane may have running biases (i.e., biases that depend 
on the source size) related to the different resolutions achieved 
in the different images (i.e., the structure is better resolved as the 
source expands). One way to avoid such biases is to use a dy- 
namic be am, i.e., a similar rati o between source size and beam 
size (e.g. iMarcaide et al]|1997h at each epoch. Here we use sev- 
eral different methods to determine the source size. The results 
for each method are summarized in Table |2] 

4.LL Concentric Rings 

First, we have measured the intensity in concentric rings (of 0.05 
mas width) using the AIPS task IRJNG. Examples of intensity 
profiles for both observations are shown in Fig. [3] together with 
fits to the data. One can define the radius Rn% where the inten- 
sity falls to n% of the peak intensity, with an arbitrary value for 
n. Since Rn% can be multi-valued (i.e., the inner and the outer 
radius), we always use the outer radius. Fig. [3] shows also the 
Rioo% (i-e. the peak itself) and Rso*, values. A value of n-5Q% 
has the advantage that the profile is much steeper and the posi- 
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Fig. 3. Intensity of concentric rings of widths 0.05 mas in the 
images from 2008 May 03 and 2009 April 08. The flux densities 
in the second observation were scaled by a factor of 5. The ver- 
tical lines mark the position of the peak and the radius where the 
intensity is 50% of the peak level. 



tion can be better determined than at the peak (where the profile 
is flat). In order to avoid a running bias, it is important to use in 
each epoch a similar ratio of source size to beam size (i.e. a dy- 
namic beam). This method will give a size which is equal to the 
real size of the expanding shell times an unknown factor This 
factor will be smaller for large values of n (the measured sizes 
are smaller for large values of «). Figure |4] shows the images of 
SN 2008iz at both epochs and an image of the second epoch with 
a dynamic beam together with the two R5o% source sizes. 

It is important to note that the intensity profile is sensitive to 
the position of the central pixel of the concentric rings. A posi- 
tion offset would smear out the intensity profile. This can be seen 
in Fig. |5] where we plot the width of the intensity profile for dif- 
ferent position offsets in the first observation. As expected, the 
width shows a strong dependence on the position offset, and the 
position of the center of the shell can be estimated. The position 
of the center is located in the first epoch at 09*'55™5P.55026, 
+69°40'45".7913 (J2000). The uncertainty in the absolute posi- 
tion is dominated by the uncertainty of the position of the phase 
referencing source M81* (~0.5 mas). 



4.L2. Common Point IVIetliod 

A different method that uses the concept of a dynamic 
beam in a natural wa y is the Common-Point Method (CPM, 
IMarcaide et al.ll2009ah . This method relies on the existence of 
a point in the radial profile of the supernova structure that re- 
mains unaltered under small changes in the convolving beam. 
This point is closely related to the source size (indeed, the ratio 
between the radial position of this point and the s ource radius 
is ~ 1; see Appendix A of IMarcaide et alj|2009al for details). 
In Fig. |6] we show the radial intensity profiles of SN2008iz, 
computed from the azimuthal averages of the CLEAN models 
convolved using different beams. The "common points" in the 
profiles can be clearly seen in the figure. 
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Fig. 4. HSA images at 22 GHz of SN2008iz in May 2008 (top) 
and April 2009 (middle and bottom). The bottom image is con- 
volved with a dynamic beam of 0.99 mas (i.e. the ratios of beam 
size and source size are identical in the top and bottom images). 
The rings denote the best fitted positions and radii from CPM in 
both epochs, i.e. 320 and 1580 //as. 

4.1.3. Model Fitting 

Finally, we estimated the size from model-fitting to the visi- 
bilities, using a simplified model of the supernova radio struc- 
ture. This method may have undetermined running biases if the 
model fitted to the visibilities is not a good representation of the 
true supernova emission structure. We fitted the visibilities us- 
ing a model of a spherical shell of 30% fractional width. This 
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Fig. 5. The width of the intensity profile as a function of position 
offset (relative to a reference pixel) in the observation on 2008 
May 03 One pixel corresponds to 0.05 mas. Shown are position 
offsets in right ascension (-I-) and declination (x). 
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Fig. 6. Solid lines: radial intensity profiles of SN2008iz at the 
two VLBI epochs, computed using different convolving beams 
(1.0, 1.1, 1.2, 1.3, and 1.4timesthe source radius at each epoch). 
The radial positions of the common points at each epoch (see 
Section |4. 1.2b are marked with dashed lines. 



is the model that best describes the structure of SN 1993 J in 
iMarcaide et al.l (l2009ah . A smaller frac tional width between 20 
and 25%, as prefered by other authors ("B artel et al."2000l l2002l: 
Bietenholz et al. 2003) or predicted by Chevalier ( 1982), could 
lead to different source sizes. However, the current data, based 
only on two epochs at one frequency, do not allow a relaible 
statement about the true fractional width of SN2008iz. For the 
model-fitting, we used a modified version of the subroutine mod- 
ELFiT in the program difmap (Shepherd et al. 1995), which fits a 
parameterized shell model to the real and imaginary parts of the 
visibilities. The use of the real and imaginary parts of the visi- 
bilities (instead of amplitudes and phases) is more robust from a 
statistical point of view, especially in the case of a low signal-to- 
noise ratio, since the noise involved in the data is gaussian-like 
(as assumed in the modelling algorithms based on a mini- 
mization). 
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Table 2. Source sizes, expansion speeds, a lower limit for the 
expansion index m, and explosion dates (in 2008) assuming m-l 
and m-0.E9 for the different methods outlined in Section|4T| 



Method 


R2OO8 


R20051 




to 


m 


to 




[juas] 


[jtias] 


km s"' 


(m=l) 




(m=0.89) 


R|00% 


149 


1008 


15800 


Mar 05 


>0.86 


Mar 19 


R50% 


307 


1462 


21200 


Feb 02 


>0.70 


Feb 22 


R25% 


380 


1676 


23800 


Jan 24 


>0.66 


Feb 14 


CPM 


320 


1580 


23100 


Feb 07 


>0.71 


Feb 25 


MODELFIT 


320 


1390 


19600 


Jan 22 


>0.65 


Feb 12 



4.2. Expansion speed 

Using the R5()% value as physical size, we get a radius of 307 
//as on 2008 May 03 and 1462 fias on 2009 April 08. The mea- 
sured sizes correspond to an average expansion speed of 1243 
fias yr"' or Vexp=21200 km s"'. Extrapolating this expansion 
back to a radius of zero and assuming constant expansion (r oc 
t'" with an expansion index m-l) leads to an explosion date to 
of 2008 February 02. However, a significant deceleration (i.e. 
m < 1) would shift the explosion date to a later time. Since the 
supernova was first seen in a VLA observation on 2008 March 
24, one can use this to give a low er limit for the expansion in- 
dex m of 0.7. lMarchih et alJ(l20ld ) model the 5 GHz light curve 
of SN 2008iz and find evidence of a modest deceleration with an 
expansion index m=0.89 and an explosion date of 2008 February 
18 (±6 days). Extrapolating the expansion seen on the VLBI im- 
ages backward using this expansion index, yields an explosion 
date of 2008 February 22. This value for the explosion date is 
in reasonable agreement with the estimate from the 5 GHz light 
curve. 

The resulting source sizes for the different methods are sum- 
marized in Table[2] Most methods give similar expansion speeds 
between 19600 and 23800 km s"'. Since the Rioo% values give 
a significantly lower expansion speed and a very late explosion 
date, we conclude that the Rioo% values underestimate the true 
source sizes significantly. Hence, we do not consider these val- 
ues for our further analysis. Based on the information of the two 
VLBI images, we thus conclude that the explosion occured be- 
tween 2008 January 22 and 2008 March 24. We get a lower limit 
for the expansion index of m >0.65, and an average expansion 
speed in the range of 19600 - 23800 km s"' . 

Since the CPM source sizes should be very close to the real 
sizes, the expansion speed is very close to the average value of 
all methods (~22000 km s"', excluding Rioo% which is clearly 
an outlier), and the explosion date is in very good agreement 
with the light curve modeling, we adopt these values. The formal 
errors in our various fits are exceedingly small (of the order of 
20 fias). Since the systematic errors from the different methods 
are much larger, the spread in velocities (+1650 km s"') and 
explosion dates (+ 7 days) reflects our uncertainties. 

This expansion speed is much higher than the values reported 
in lBrunthaler et al] (l2009bl) . This discrepancy has three reasons: 
i) The lower value was based on a preliminary data reduction, 
without the data from the phased VLA and without the geodetic 
block corrections, ii) The radius of the brightest emission was 
used as an approximation for the source size. However, as shown 
here, this method underestimates the true expansion speed, iii) 
The sizes were measured by hand from the images and this gave 
an overestimated size in the first epoch, where the source was 
still very compact. 
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Fig. 7. Cross-correlation coefficient, C(a, 6), between the images 
of the two VLBI epochs (scaling the first epoch to the size of the 
second), as a function of the coordinates of the explosion center, 
(a, 6), taking as reference the position of the maximum C(a, 6) 
(see Sect.l43Tl. 



4.3. Expansion center, self-similarity, and anisotropic 
expansion 

The ring is not symmetric and the brightest region in both im- 
ages is the south-western part of the shell. This could be caused 
either by an asymmetry in the explosion or in the circumstellar 
medium (CSM), i.e. a clump with higher density in the south- 
west . A clumpy CSM is very common for red supergiant stars 
(e.g. lSmith et al. 2009) which are likely progenitors for this su- 
pernova. Furthermore, a comparison of the 5 GHz light curve 
and t he available 22 GHz data shows evidence for a clumpy 
CSM (" MarchiU et al] |20T^. 

As mentioned in section 14.1.11 we also estimated the cen- 
tral position of the ring. We do not detect any significant shift in 
the position, with an upper limit of ~ 50 //as in right ascension 
and 100 /(as in declination. To verify the relative astrometric ac- 
curacy of the positions, we also imaged the three extragalactic 
background sources 0945+6924, 0948+6848, and 1004+6936. 
The average position change between the two observations of 
M8 1 * relative to these quasars is 1 + 15 //as in right ascension and 
22+53 //as in declination, where the uncertainties are the stan- 
dard deviatiorQ. Thus, we can rule out that any significant posi- 
tion shift is introduced by jet motion in M81*. We note that the 
angular separation of M82 to M8 1 * is not larger than the angu- 
lar separation between M81* and the three background quasars. 



^ These values give 1 cr upper limits on the proper motion of M81 of 
260 km s"' in right ascension and 1270 km s"' in declination 
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Hence, systematic errors in the astrometry should be similar for 
all sources. 

According to lChevaheii (Il982h . the expansion of a supernova 
is self-similar, i.e., the source structure remains unaltered regard- 
less of the source size. This model of a self-si milar expansion 
has b een ext ensively tested w ith SN 1993J (e.g. iMarcaide et all 
11997 . 20Q'9a;' Bartel et al.l2002i) . Although some deviations from 
self-similarity were found in the expansion of SN 1993 J, this su- 
pernova kept its structure nearly self-similar during more than a 
decade. If the expansion of SN 2008iz was also self-similar be- 
tween our two VLBI epochs, the VLBI images obtained should 
be equal (regardless of the global flux density scale) if we scale 
the first image to the size of the second one and use the same con- 
volving beam for both images. However, the result of expanding 
the first VLBI image, to compare to the second one, depends on 
the point in the image that is chosen as the center of expansion. 
Therefore, a direct comparison between the images, to check the 
self-similarity in the expansion, is not possible if the coordinates 
of the expansion center of the supernova are not known. 

To determine the coordinates of the expansion center of the 
supernova independently from the method described above, we 
computed different expanded versions of the image of the first 
VLBI epoch using different centers of expansion. The scaled im- 
ages were obtained by scaling the positions of the CLEAN com- 
ponents of the model, according to the ratios of the supernova 
sizes reported in the previous section. We then used the same 
beam to convolve the CLEAN model of the second epoch and 
the resulting (scaled) CLEAN models of the first epoch. Finally, 
we compared the resulting image of the second epoch with the 
resulting (scaled) images of the first epoch by computing the 
cross correlation between both images using: we mean 

2, 4 X/; (a, 5) 

Cifl, 6) = ^ \ 2 (2) 

where I'^ is the ith pixel of the image of the second epoch and 
I[{a,6) is the /th pixel of the scaled version of the image of 
the first epoch, taking the point (a, S) as the expansion center 
Assuming a self-similar expansion, the coordinates of the max- 
imum value of C{a,5) are an estimate of the position of ex- 
pansion center of the supernova, which we identify as the co- 
ordinates of the explosion. In Fig. [T] we show the cross cor- 
relation of the images computed for different expansion cen- 
ters, {a, 6), expanding the image of the first epoch according 
to the size estimated with the CPM. The maximum correla- 
tion (which corresponds to our estimate of the explosion cen- 
ter) takes place at the coordinates a - 09*' 55™ 5 P. 55025 and 
6 = -1-69° 40' 45". 79133. These coordinates do not change by 
more than 10//as if we use, instead, the size estimated with 
model-fitting to expand the image of the first VLBI epoch. This 
position agrees to the position given in Sect. I4.L1I within 60 
//as and verifies that the position of the expansion center has not 
changed between the two epochs. 

Additionally, the value of C(a, 6) at the maximum is a mea- 
sure of the degree of self-similarity (and/or anisotropy) in the ex- 
pansion. The maximum cross-correlation of the images is 0.98, if 
we use the CPM size to expand the image of the first epoch, and 
0.97, if we use instead the model-fitting size. These values are 
very close to 1, which is the case of a perfect self-similar expan- 
sion. Therefore, we conclude that the expansion of SN2008iz 
was self-similar to a high degree between our two VLBI epochs. 




from lWeiBetaLl (1200 lb with the positions of SN2008iz (stai-) 
and the MERLIN transient (cross). The contours correspond to 
200, 400, 600, 800, 1200 1600, 2000, and 2400 K kms"'. The 
resolution of the observation is 1.5" x 1.4" (~25 pc). 



5. Discussion 

5.1. Column density and extinction 

The non-detection of SN2008iz in the optical, infrared, and 
X-rays indicates that it exploded inside or behind a very 
dense cloud. Inde ed, the '^CO (J=2— >1) line intensity map in 
IWeiBetal.l(l2001h shows a prominent cloud exactly at the posi- 
tion of the supernova (Fig. [8]). The line intensity at the position 
of SN 2008iz is ~ 1800 K kms"' . Using the Galactic conversion 
factor Xco = 1.6 x 10^*' cm"^ (K kms"')"', this corresponds 
to a H2 column d ensity A^(H2) of ~ 29 x 10^^ cm"^. However, 
IWeiB et al.l (1200 ih find much smaller and spatially variable con- 
version factors from radiative transfer calculations, that lead to 
smaller H2 column densities. At the position of SN2008iz, the 
conversion factor is ~ 0.3 x 10^" cm"- (K kms"')"' (their Fig. 
10). This leads to a H2 column density of ~ 5.4 x 10^^ cm"^. 
However, the CO observations were performed with a linear res- 
olution of ~25 pc. Thus it is possible, that the supernova is lo- 
cated behind a smaller cloud with much higher column density. 
The '^CO line intensity at the position of the MERLIN transient 
is almost 4 times smaller (~500 K kms"'). 

Taking the latter value (5.4 x 10^^ cm"^) for N(H2) and as- 
suming that all of the column is between SN 2008iz and us, 
we derive a visual extinction. Ay, of 24.4 mag. Here we have 
used the relation between optical extinction and hydrogen nu- 
cleus column density, N(H), derived by Giiver & Ozel (20091) 
from X-ray absorption data of a sample of Galactic SNRs and 
assumed that all the hydrogen is in molecular form, i.e., N(H2) 
= 2 N(H). Such a high value of the extinction would explain the 
lack of an optical counterpart. 

The deriv ed extinction is muc h higher than toward SN 1 993 J, 
for which Ric hmond et al ] (fT994l) discuss values of Ay = 0.7 an d 
1 .0 mag. With the extinction law given bv^ Cardelli et al.l (Il989l) . 
Ak - 0.114 Ay , we c alculate 2.8 mag for the A'-band (2.2 fim) 
extinction A^. iFraser et al.l (12009) report an even higher /T-band 
extinction of up to 1 1 mag based on their non-detection (3cr up- 
per limit on the absolute /T-band magnitude of -5 mag) and the 
assumption that the infrared light curve of SN 2008iz behaves 
similar to the one of SN 1993J. Thus, SN2008iz is either very 
weak in the infrared compared to SN 1993 J, or behind a smaller 
but even denser cloud than estimated here. 
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Table 3. Comparison between SN 1993 J and SN 2008iz. 



Property 




SN1993J 


SN2008iz 


spectral index a 




-0.99" 


-1.08" 


LsGHz 


[lO^'ergs-' Hz-'] 


1.5'' 


2.5'' 


tpeak-to 


[days] 


180' 


-120'' 


VVLBI 


[kms-'] 


14900" 


21200* 


Lx-rayatt~220days 


[1038 ergs-'] 


~8-f 


< 15'' 



Refere nces: a) van Dyk et al. (1994); b) this work; c) Weileretal. 
(Il998h; dl iMarchili et al. (2010); e) iMarcaide et all 11995.) ; f) 
IZimmermann & Aschenbactj ( l2003h . 



5.2. Expansion velocity and syncinrotron self-absorption 

IChevalie3 ( Il998h proposed a way to estimate the mean expan- 
sion velocity of a supernova based on the radio-light curve and 
assuming that synchrotron self-absorption (SSA) dominates the 
light curve at the peak. If SSA is not the dominant absorption 
process at the peak of the li ght curve, then t he estimate of the 
expansion velocity from the ChevalieJ (Il998l) model is a lower 
bound to t he real expansion velocity of the radio shell. Following 
IChevalieJ (199 8). the estimated mean expansion velocity at the 
peak of the radio-light curve, assuming dominant SSA, is: 



yssA(kms-i) = 5.3786 x 10^ (J34>r''"'' ( x 



D \ 



IMpc 



^Q.95 ^ 

1 1 1.0 GHz/ 



Iday 



1 



(3) 



where /3 is the ratio of the relativistic particle energy density to 
the magnetic field energy density (if we assume ene rgy equipar- 
tition,y8 = i(i+k) ' where k ranges from 1 to 2000, see iPacholczvM 
\T9Td . chapter 7), is the filling factor of the emitting region to 
a sphere (we assume a shell of width equal to 30% of the outer 
radius, which yields <p = 0.66), Fp is the flux density at the peak, 
D is the distance, v is the observing frequency, tp is the super- 
nova age at the pe ak, and w e use a spectral index a = -1 (see 
Eqs. 11 and 13 of Chevalier 1998). 

For SN2008iz, Eq. [3] yields a mean expansion velocity (de- 
pending on k) in the range 8 100- 1 1800km s ' at the peak of the 
5 GHz radio light curve. Using an expansion index of 0.89, these 
velocities translate into a mean expansion velocity in the range 
7000 - 10300km S-' at the epoch of 27 April 2009. This range 
of velocities is a factor ~ 2 smaller than the velocities estimated 
from our VLBI observations, thus indicating that, in contrast to 
the case of SN1993J, SSA may not be an important absorption 
mechanism in the SN200 8iz radio emis s ion. This is also con- 
sistent with the results in llMarchili et al.l (I201Q) who were able 
to model the radio light curve of SN2008iz assuming that SSA 
effects are much smaller than free-free absorption (FFA) during 
the whole supernova expansion. 

5.3. Comparison with other type II radio supernovae 

Estimates of the expansion velocities of other type II ra- 
dio supernovae have been estimated from VLBI observations, 
and very different results have been obtained. For instance, 
the mean expansion velocity of SN 1979C during the first 
year after explosion is estimated to be ~ 10000 - 11000 
km S-' (Bartel & Bietenholz 2003; Marcaide et al. 20P91l); for 
SN 1986J, a velocity of ~ 1470 km s^' was obtaine d by 
iPerez-Torres et all (l2002bl) . while iBietenholz eTaP (l2002h find 



20000 km s ' 3 month after the explosion: IStavelev-Smith et al.l 
(1 19931) report a mean expansion speed of ~ 35000 km s ' dur- 
ing the first years for SN 1987 A before it slowed down to ~ 
4800 km s '; f or SN2004et, the expa nsion velocity was > 
15700 km s ' dlVIartf-Vidal et al.1 l2007l) : and for SN2008ax, 
an expansion velocity as large as 52000 km s ' was obtained 
(iMartf-Vidal et al.l l2009l) . Estimates of the expansion veloci- 
ties of other supernova remnants in M82 (the host galaxy of 
SN2008iz), have been al so reported, which r ange between ~ 
1500 and 11000 km s^^ dBeswick et al.ll2006h . These later ve- 
locities are much higher (a factor of 3 - 22) than the pr edicted 
velocities from the model of^ Chevalier & Fransson (2001), based 
on the high pressure expected in the interstellar medium (ISM) 
of M82. The expansion velocity reported in this paper for 
SN2008iz is indeed a facto r ~ 4 large r than the predicted ve- 
locities in Chevalier & Fran ssonI (1200 ih . although of the same 
order of magnitude than the velocities reported in lBeswick et alj 
(2006.) for the other remnants in M82, and the typical velocities 
of t he other type II sup ernovae observed to date. 

IWeileretan(ll998h find a correlation between peak radio lu- 
minosity at 5 GHz and the time between the explosion and the 
peak in the 5 GHz light curve for type II supernova e. The 5 GHz 
light curve of SN 2008iz from lMarchili et all (1201 Oh gives a peak 
luminosity of ~ 2.5 x 10^^ erg s"' Hz^' at ~120 days after the 
explosion. These values are well within the scatter of the corre- 
lation. Thus it seems plausible that SN 2008iz is also a type II 
supernova. 

Since SN2008iz and SN 1993J are located at very similar 
distances, this allows a detailed comparison between these two 
supernovae. Several properties of both supernovae are summa- 
rized in Table |3] The radio spectral indices, the peak radio lu- 
minosities, rise times, and early VLBI expansion velocities are 
similar (considering that the rise times and peak radio luminosi- 
ties can vary by several orders of magnitudes for type II radio 
supernovae). The non-detection in X-rays of SN2008iz can be 
attributed to absorption by the dense molecular cloud seen in the 
CO data. 



6. Summary 

In this paper we presented the first VLBI images, a VLA radio 
spectrum from 1 .4 to 43 GHz, and Chandra X-ray observations 
of SN 2008iz. Our main results are: 

- The VLBI images, separated by ~ 1 1 month show a shell- 
like structure expanding with a velocity of ~ 23000 km s"'. 

- The inferred expansion speed is a factor of 2 higher than 
expected if SSA dominates the light curve. This indicates 
that SSA is not important for the radio emission. 

- The most likely explosion date is in mid February 2008, but 
not earlier than January 22 and not later than March 24. 

- We find no evidence for an asymmetric explosion, but a 
high degree of self-similarity between the two VLBI obser- 
vations. 

- The VLA radio spectrum is well fitted by a broken power 
law with a turnover frequency of 1 .5+0. 1 GHz, and a spectral 
index of -1.08 ± 0.08 in the optically thin part. 

- SN2008iz is located behind (or inside) a large molecular 
cloud with a H2 column density of 5.4x 10^^ cm"^ (on a scale 
of 25 pc), corresponding to a visual extinction. Ay, of 24.4 
mag. 

- Due to the high column density, we obtain only upper limits 
on the X-ray luminosity of 1.5x10'*' and 1.5x10-''' erg s"' 
~75 and ~200 days after the explosion, which is consistent 
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with the X-ray luminosity of SN 1993J at similar ages. We 
also obtain an upper limit for the X-ray luminosity of the 
second radio transient of 1.2x10^^ erg s"^ (in 0.3-2.4 keV). 
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